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A coupled wave analysis is given of the Bragg diffraction of light by 
thick hologram gratings, which is analogous to Phariseau's treatment of 
acoustic gratings and to the "dynamical" theory of X-ray diffraction. The 
theory remains valid for large diffraction efficiencies where the incident 
wave is strongly depleted. It is applied to transmission holograms and to 
reflection holograms. Spatial modulations of both the refractive index and 
the absorption constant are allowed for. The effects of loss in the grating and 
of slanted fringes are also considered. Algebraic formulas and their nu- 
merical evaluations are given for the diffraction efficiencies and the angular 
and wavelength sensitivities of the various hologram types. 

I. INTRODUCTION 

Holographic recording in thick media ("volume recording") is of 
particular interest for high-capacity information storage, 1-3 for color 
holography 4 and for efficient white-light display of holograms. 5-9 The 
high efficiency of light conversion which is attainable with thick di- 
electric holograms is also important for microimaging, and it may make 
it practical to use holographic optical components (for example, gratings 
or fly's eye lenses) in a variety of optical systems. 

In thick holograms it is light diffraction at or near the Bragg angle 
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which leads to efficient wavefront reconstruction. This is true for both 
transmission and reflection holograms, and both types are considered 
in this paper. The (volume) record of the holographic interference 
pattern (fringe pattern) usually takes the form of a spatial modulation 
of the absorption constant or the refractive index of the medium, or 
both. Modulations of the absorption constant are produced in con- 
ventional photographic emulsions and in photochromies, while newer 
materials, like dichromated gelatin 10,11 lithium niobate, 12 or photo- 
polymer materials 13 yield modulations of the refractive index. 

This paper considers the properties of all these types of thick (or 
"deep") holograms. Of particular interest is their efficiency of convert- 
ing light into the useful reconstructed wave (diffraction efficiency) and 
the angular dependence of this diffraction efficiency as the incident 
light deviates from the Bragg angle. We are also interested in the wave- 
length dependence and in the way the diffraction properties are changed 
in the presence of loss or a slant of the fringe pattern with respect to 
the surface of the recording medium. 

Leith and his associates, and Gabor and Stroke have already con- 
sidered some of the properties of thick holograms, in particular the 
angular and the wavelength dependence of the diffracted light. 14,15 Their 
theories are essentially linear or perturbational theories which use the 
Kirchhoff integral or the first Born approximation with the basic 
assumption that the incident light wave is not disturbed by the dif- 
fraction process. Their results are valid as long as this assumption is 
good. For high diffraction efficiencies (like 90 percent) the incident 
wave is strongly depleted and another approach is called for. One such 
approach is to use electronic computers to solve the relevant compli- 
cated electromagnetic problem accurately. Results of such computations 
are available for special cases. Klein, Tipnis, and Hiedemann have com- 
puted data for light diffraction by ultrasonic waves, 1G ' 17 and Burckhardt 
has reported results for dielectric hologram gratings. 18,19 The method 
of Bathia and Noble 20 is another approach in which they employed 
integral equations to analyze acoustic diffraction of light. 

Yet another approach is the use of a coupled wave theory, which is 
the subject of this paper. Such a theory can predict the maximum 
possible efficiencies of the various hologram types (results which one 
cannot hope to obtain from linear theories), and the angular and wave- 
length dependence at high diffraction efficiencies. Following Phariseau, 21 
coupled wave theories have been successfully used in the treatment of 
light diffraction by acoustic waves 22 and by electrooptic gratings 23 
where very similar diffraction processes are at work as in holography. 
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Closely related to the diffraction in thick holograms are also the dif- 
fraction of electrons in lattices and the diffraction of X-rays in crystals. 
The dynamical theory of X-ray diffraction 24 is also a theory of coupled 
waves and its application to holography has already been suggested. 25 

We have earlier reported some of the results and an outline of the 
coupled wave theory for hologram gratings. 26,27 Here we propose to 
give further results and a more detailed account of the basic assumptions 
and the analysis. We give analytic formulas for the various hologram 
types as well as numerical evaluations which include results on the 
angular sensitivities and the influence of loss and slant. 

For simplicity the analysis is restricted to the holographic record 
of sinusoidal fringe patterns which we call hologram gratings. To some 
degree a more complicated hologram can be regarded as a multiplicity 
of such hologram gratings. 

II. COUPLED WAVE ANALYSIS 

2.1 Derivation of the Coupled Wave Equations 

The coupled wave theory assumes monochromatic light incident on 
the hologram grating at or near the Bragg angle and polarized per- 
pendicular to the plane of incidence.* Only two significant light waves 
are assumed to be present in the grating: the incoming "reference" 
wave R and the outgoing "signal" wave S. Only these two waves 
obey the Bragg condition at least approximately, the other diffraction 
orders violate the Bragg condition strongly and are neglected. They 
should be of little influence on the energy interchange between S 
and R. The last assumption limits the validity of the coupled wave 
theory to thick hologram gratings. Section 6 gives a more detailed 
discussion of this limitation. 

Figure 1 shows the model of a hologram grating which is used for 
our analysis. The z-axis is chosen perpendicular to the surfaces of the 
medium, the .r-axis in the plane of incidence and parallel to the medium 
boundaries and the y-axis perpendicular to the paper. The fringe 
planes are oriented perpendicular to the plane of incidence and slanted 
with respect to the medium boundaries at an angle <f>. The fringes 
are shown dotted. The grating vector K is oriented perpendicular to 
the fringe planes and is of length K = 2tt/A, where A is the period 
of the grating. The same average dielectric constant is assumed for 
the region inside and outside the grating boundaries. The angle of 
incidence measured in the medium is 0. 

* A generalization to parallel polarization is given in the appendix. 



2912 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1969 




|* d -*| 

Fig. 1 — Model of a thick hologram grating with slanted fringes. The spatial 
modulation of n or a is indicated by the dotted pattern. The grating parameters are: 
— angle of incidence in the medium, K — grating vector (perpendicular to the fringe 
planes), A — grating period, 4> — slant angle, and d — grating thickness. 

Wave propagation in the grating is described by the scalar wave 
equation 

S7 2 E + k 2 E = 0, (1) 

where E(x, z) is the complex amplitude of the y-component of the 
electric field, which is assumed to be independent of y and to oscillate 
with an angular frequency w. The propagation constant k(x, z) is 
spatially modulated and related to the relative dielectric constant 
e(x, z) and the conductivity <j{x, z) of the medium by 



? 2 w 

K = — € — JCOMC 

c 



(2) 



where c is the light velocity in free space and n is the permeability 
of the medium which we assume to be equal to that of free space. 
In our model the constants of the medium are independent of y. The 
fringes of the hologram grating are represented by a spatial modulation 
of « or a: 
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e = e + e, COS (K-x) 
a = o-„ + <T, COS (K-x) 



(3) 



where «i and <j y are the amplitudes of the spatial modulation, e is 
the average dielectric constant and <r the average conductivity, e and 
o- are assumed to be modulated in phase. To simplify the notation 
we have used the radius vector x and the grating vector K 
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x = 



Equations (2) and (3) can be combined in the form 

A; 2 = /3 2 - 2ja(3 + 2 K p(e iKx + e~ mx ) 

where we have introduced the average propagation constant /3 and 
the average absorption constant a: 

& = 2ir(« )*A; « = jwar^Ceo) 1 , (5) 

and the coupling constant k was defined as 



1 2 



K = 4 



6i/(eo) J - i/*c«ri/(e o )0 • 



(6) 



This coupling constant describes the coupling between the reference 
wave R and the signal wave S. It is the central parameter in the coupled 
wave theory. For k = there is no coupling between R and S and, 
therefore, there is no diffraction. 

Optical media are usually characterized by their refractive index 
and their absorption constant. We also find it convenient to use these 
parameters if the following conditions are met 



2tt7i/\ » a; 2mi/X » «i , n » rij , 



(7) 



which is true in almost every practical case. Here n is the average 
refractive index, and n x and a x are the amplitudes of the spatial modula- 
tion of the refractive index and the absorption constant, respectively 
[compare equation (3)]. X is the wavelength in free space. Under the 
above conditions we can write with good accuracy 



= 2irn/\ 



(8) 
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and for the coupling constant 



K = TTH-j/X — j'ai/2. 



(9) 



The spatial modulation indicated by ?ii or a { forms a grating which 
couples the two waves R and S and leads to an exchange of energy 
between them. We describe these waves by complex amplitudes R (z) 
and S(z) which vary along z as a result of this energy interchange or 
because of an energy loss from absorption. The total electric field in 
the grating is the superposition of the two waves: 

E = R(z)e- iex + S(z)e-" x . (10) 

The propagation vectors 9 and d contain the information about the 
propagation constants and the directions of propagation of R and S. 
9 is assumed to be equal to the propagation vector of the free reference 
wave in the absence of coupling, d is forced by the grating and related 
to 9 and the grating vector by 



d = p - K 



(11) 



which has the appearance of a conservation of momentum equation. 
9 and d have been chosen to conform as closely as possible with our 
picture of the physical process of the diffraction in the grating. If the 
actual phase velocities differ somewhat from the assumed values, then 
these differences will appear in the complex amplitudes R (z) and S(z) 
as a result of the theory. 

Figure 2 shows the vectors of interest and their orientation. The 
components of 9 are p x and p z which are given by 







sin | 





= 





-P*. 




1 cos 0, 



9 = 



From this and equation (11) follow the d-components <r x and a t 

K 



(12) 



d = 





10" .-J 



= jS 



sin 6 — sin <f> 


cos 6 — — cos <f> 



(13) 



The vector relation (11) is shown in Fig. 3 together with a circle 
of radius /3. The general case is shown in Fig. 3a, where the Bragg 
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Fig. 2 — o and d, the propagation vectors of the reference wave R and the signal 
wave S, and their relation to the grating vector K. The obliquity factors c R and c s 
are indicated. 

condition is not met and the length of d differs from /3. Figure 3b shows 
the same diagram for incidence at the Bragg angle O . In this special 
case the lengths of both, p and d are equal to the free propagation 
constant /3, and the Bragg condition 



cos (0 - 0) = K/20 



(14) 



is obeyed. 

For a fixed wavelength the Bragg condition is violated by angular 





Fig. 3 — Vector diagram (conservation of momentum) for (a) near and (b) exact 
Bragg incidence, 
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deviations A from the Bragg angle O ■ For a fixed angle of incidence a 
similar violation takes place for changes AX from the correct wave- 
length X • We write 

e = d + A0, (15) 

and 

X = X + AX, 

and assume in the following that the deviations A0 and AX are small. 
Angular changes A0 have very similar effects on the behavior of the 
grating as wavelength changes AX, and there is a close relation between 
the angular sensitivity and the wavelength sensitivity of thick hologram 
gratings. We get an idea of this relationship by differentiating the 
Bragg condition (14), from which results 

^ = £/4imsinfo- 0„). (16) 

aX 

The — X connection shows up in the dephasing measure t? which 
appears in the coupled wave equations and which is defined by 

d ee OS 8 - cr 2 )/2/3 = K cos (0 - 0) - ¥- X (17) 

and which has been expressed in this form using equation (13). A Taylor 
series expansion of equation (17) yields the following expression for # 
which is correct to the first order in the deviations A0 and AX: 

d- = Ad-K sin (0 - O ) - A\-K 2 /4m. (18) 

Note that the deviations A0 and AX which produce equal dephasing 
# are related by equation (16). 

We are now ready to derive the coupled wave equations. We combine 
equations (1) and (4), and insert the expressions of (10) and (11). 
Then we compare the terms with equal exponentials (e~' g " and e~ ia,t ) 
and arrive at 

R" - 2jR' P , - 2jc0R + 2k0S = (19) 

and 

S" - 2jS'<r z - 2ja$S + (/3 2 - a 2 )S + 2 K pR - 0, (20) 

where the primes indicate differentiation with respect to z. The waves 
generated in the directions of p + K and d — K are neglected, together 
with all other higher diffraction orders. In addition we assume that the 
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energy interchange between S and R is slow and that energy is absorbed 
slowly, if at all. This allows us to neglect R" and S". We will check 
the results of the theory later for a more detailed justification of this 
last step. We can now introduce equation (18) and rewrite the above 
equations in the form 

c R R' + aR = -jkS (21) 

c s S' + (a + jd)S = -JkR. (22) 

These are the coupled wave equations which are the basis for our 
analysis. The abbreviations c R and c s stand for the expressions 

c R = p,//3 = cos e , 23 s 

Cs = Gz/fi — cos — — COS <f>. 

Our physical picture of the diffraction process is reflected in the coupled 
wave equations. A wave changes in amplitude along z because of coupling 
to the other wave (kR, kS) or absorption (aR, aS). For deviations from 
the Bragg condition S is forced out of synchronism with R and the 
interaction decreases (&S). 

The energy balance of the coupled-wave model is described by the 
relation 

(c R RR* + c s SS*Y + 2a(RR* + SS*) + j(k - k*)(RS* + R*S) = 

(24) 

where the asterisk denotes a complex conjugate. This is easily derived 
from equations (21) and (22) by multiplying them with R* and S*, 
respectively, and adding the results together with the complex con- 
jugate results. The presence of the obliquity factors c R and c s in the 
first part of equation (24) indicates that it is the power flow of the 
two waves in the z direction that enters the energy balance. In the 
absence of ohmic loss this power flow is conserved. The second and the 
third part in the equation describe the energy loss resulting from ab- 
sorption in the grating. They correspond to the relevant terms of a-EE*. 

2.2 Solution of the Coupled Wave Equations 

It is straight forward to obtain the general solution of the coupled 
wave equations, which is 

R{z) = r x exp (txz) + r 2 exp (y 2 z) (25) 

S(z) = Si exp (7,2) + s 2 exp (y 2 z) (26) 
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where the r,- and s,- are constants which depend on the boundary condi- 
tions. To determine the constants 7,- we insert equations (25) and (26) 
into the coupled wave equations and obtain 

(c«T. + «>, = -JKSi (27) 

i = 1,2 

(fisyi + a + j&) Si = -JKr< . (28) 

After multiplying the equations with each other we get a quadratic 
equation for t< 

(c R7i + a)(c syi + a + m = -k\ (29) 

with the solution 

_1 /«. , a, •*_ 

2 \Cr Cg Cg. 

±ir(----i-T-4— T- (30) 

At this point we divert briefly from the main derivation, because 
now we have the means to check the validity of neglecting R" and 
S" in Section 2.1. This step is justified if the conditions R" « £,R', 
and #" <& <r t S' are obeyed. In view of equations (25) and (26) this 
will happen if y,- <3C /S. According to equation (30) the above requirement 
is met if Ad <5C 1 and if the inequalities of equation (7) are satisfied 
(which is usually the case). 

Continuing the coupled wave analysis, we have to determine the 
constants r,- and s, . To do this we have to introduce boundary condi- 
tions into our model. These are different for transmission holograms 
and for reflection holograms. Figure 4 gives an indication of this. For 
both hologram types the reference wave R is assumed to start with 
unit amplitude at z = 0. It decays as it propagates to the right and 
couples energy into #. In transmission holograms the signal S starts 
out with zero amplitude at z = and propagates to the right (c s > 0). 
In reflection holograms the signal travels to the left (c s < 0) and it 
starts with zero amplitude at z = d. 

Let us first analyze transmission holograms where c s > 0. Here, 
the boundary conditions are 

R(Q) = 1, S(0) = (31) 

as discussed before. If we insert these boundary conditions into equa- 
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Fig. 4 — Wave propagation in (a) transmission and (b) reflection holograms. 
The reference wave R decays while it propagates to the right. In (a) the signal S 
travels to the right and gains with z, while in (b) S travels to the left and gains with 
decreasing z. The shading indicates the orientation of the fringes. 

tions (25) and (26), it follows immediately that 

r, + r 2 = 1, 
and (32) 

Si + s 2 = 0. 

Combining these relations with equation (28) we obtain 

Si = — s a = -jn/cshx — 72). (33) 

Introducing these constants in equation (26) we arrive at an expres- 
sion for the amplitude of the signal wave at the output of the grating 



S(d) m j 



Cs (7i - 7 2 ) 



(exp (y 2 d) - exp (7,*/)). 



(34) 



This is a general expression, which is valid for all types of thick trans- 
mission holograms including the cases of off-Bragg incidence, of lossy 
gratings and of slanted fringe planes. 

The analysis of reflection holograms follows a pattern similar to the 
above. We have c s < and boundary conditions given by 



ft(0) = 1, S(d) = 0. 



(35) 
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The output plane for the signal wave is, now, at z = 0, and S(0) is 
the output amplitude of interest. Inserting the boundary conditions in 
equations (25) and (26) yields 

ri + r 2 = 1 

and (36) 

Si exp (yxd) -f- s 2 exp (y 2 d) = 0. 

To proceed with our derivation we rewrite the above relation for s, 
and s 2 in the form 

Si(exp (y 2 d) — exp (y t dj) = (sj + s 2 ) exp (y 2 d) 

(37) 
s 2 (exp (y 2 d) — exp (y^)) = — (s, + s 2 ) exp (y x d) . 

Then we sum equation (28) for i = 1 and * =» 2 and obtain the relation 

-Mr, + r a ) = -J/c = (si + s a )(a + i*) + c«(yi«i + 7 2 s 2 ). (38) 

Using the relations (37) to substitute the sum (s x + s 2 ) for the terms 
s x and s 2 in this equation we finally arrive at the result for the amplitude 
S(0) of the output signal of a reflection hologram 



// , ■« , „ Yi exp (y 2 d) - 72 exp (y^d JX 
7 I" + ^ + Cs exp ( 7a d) - exp ( 7l d) J" 



5(0) = s, + s 2 = -i« 

(39) 

This is, again, a formula of quite general validity, including off-Bragg 
incidence, loss, and slant. 

In the following sections we discuss the behavior of transmission and 
reflection holograms in greater detail, using the general formulas de- 
rived above. In these discussions a parameter of prime interest is the 
diffraction efficiency rj, which is denned as 

r, = ^ SS* (40) 

where S is the (complex) amplitude of the output signal for a reference 
wave R incident with unit amplitude, rj is the fraction of the incident 
light power which is diffracted into the signal wave. S is equal to S(d) 
for transmission holograms and equal to S(0) for reflection holograms 
in the notation of this section. But for reasons of simplicity we omit 
the arguments in the following sections. The obliquity factors c R and 
c s appear in the above definition for the same reason they have ap- 
peared in the energy balance of equation (24): in the absence of loss 
it is the power flow in the z direction which is conserved. 
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For slanted gratings another important parameter is the slant factor c 
which is denned as the ratio between the obliquity factors 

c = c R /cs = —cos 0/cos (6 — 2<f>) 

which we have expressed here, for Bragg incidence, in terms of the 
angle of incidence 6 and the slant angle <t>. Figure 5 indicates lines of 
constant c as a function of O and <f>. For transmission holograms c is 
positive (c > 0), and for reflection holograms c is negative (c < 0). 
In the diagram transmission and reflection holograms are separated by 
the line for c = °° . 

III. TRANSMISSION HOLOGRAMS 

In this section we discuss transmission holograms in greater detail. 
We give algebraic formulas and their numerical evaluations for the 
diffraction efficiencies and the angular and wavelength sensitivities of 
dielectric and of absorption gratings. This includes results on the 
influence of loss and slant. 
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Fig. 5 — The slant factor c as a function of the angle of incidence 0o and the slant 
angle <t>. c is positive for transmission holograms and negative for reflection holograms. 
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It is convenient to write the various diffraction formulas in terms 
of parameters v and £, which are redefined for each grating type. In 
these parameters are lumped together the constants of the medium 
(to, a, TOi , ai , k), the obliquity factors (c R , c s ), the wavelength, the 
grating thickness d, and the dephasing measure #. By using v and £, 
various trade-offs become immediately apparent. 

We recall that, for transmission holograms, c s is positive and the 
output signal appears at z = d. Combining equations (30) and (34) 
we obtain a general formula for the signal amplitude S of a transmis- 
sion grating 

S = ~i(^)*-exp (-ad/c R ).e ( -sin [„» - fl»/[l - f/fl*, 
v = Kd/iCnCs)' 1 , (41) 



• -*-(£-£-'*)• 



where k is the coupling constant given in equation (9), t? the dephasing 
measure of equation (18), c R and c s are the obliquity factors of equa- 
tion (23), a is the absorption constant and d the grating thickness. 
In the above form the parameters v and £ are, in general, of complex 
value. 

3.1 Lossless Dielectric Gratings 

For completeness we give the formulas for the lossless dielectric 
grating. For the unslanted case of this grating these formulas have 
been previously obtained by several workers whose prime interest 
was light diffraction by acoustic waves. 20,21 ,28 For this grating type it 
is easy to include the effect of slanted fringes.* For the lossless dielectric 
grating we have a coupling constant k = irrij/X and a — ai — 0. Equa- 
tion (41) can be rewritten in the form 

S = -;(^)V'*sin (/ + fl»/(l + fW, 

v = Trn i d/\ic R c s ) li , (42) 

£ = &d/2c s 

where v and £ have been redefined and are real-valued. The associated 
formula for the diffraction efficiency is 



* Slant was also included in the treatment of dielectric gratings in Ref. 29. 
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r] = sin 



{v 2 + rtVd + *7A 



(43) 



For significant deviations from the Bragg condition the parameters 
v and £ are of equal order of magnitude, and we can take v as independent 
of Ad or AX without causing an appreciable change in the predictions 
of equation (43) . In this equation the angular and wavelength deviations 
are represented by the parameter £ which can be written in the form 



£ = Ad-Kds'm (<f> - 6 )/2c s 
= -A\-K 2 d/Simc s 



(44) 



by using equation (18). 

The angular and wavelength sensitivities of lossless dielectric gratings 
are shown in Fig. 6, where the efficiencies as given by equation (43) 
are plotted (normalized) as a function of £ for three values of v. The 
figure shows the sensitivity of gratings with v = 7r/4 and a peak diffrac- 
tion efficiency of tj = 0.5, with v = ir/2 and a peak efficiency of 7j = 1, 
and with v = 3ir/4 and tj = 0.5. We notice that the half-power points 
are reached for values near £ = 1.5. There is some narrowing in the 
sensitivity curves for increasing values of v, and a marked increase in 
the side lobe intensity. 




Fig. 6 — Transmission holograms — the angular and wavelength sensitivity of 
lossless dielectric gratings with the normalized efficiencies tj/ijo as a function of £. 
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Fig. 7 — Transmission holograms — the angular sensitivity of a lossy dielectric 
grating with v = ir/2 and D = 2 compared with that of a lossless dielectric grating 
(Do = 0), for O = 30° and /3d = 50. 

The above formulas include the influence of slant through the obliquity 
factors c R and c s . If there is no slant (0 = 7r/2) and if the Bragg condi- 
tion is obeyed then c R = c s = cos O and equation (43) becomes the 
well known 20 ' 21 ' 30 



t) = sin 2 {vnxd/\ cos O ). 



(45) 



By inserting the above half-power values for £ into equation (44) we 
obtain simple rules of thumb for the angular and spectral half-power 
bandwidths of unslanted gratings: 2A0j ^ A/d, 2AXj/X £3 cot 9-A/d. 

3.2 Lossy Dielectric Gratings 

Let us first study the influence of loss on the angular sensitivity 
of a dielectric grating. We assume that there is no slant (<f> = ir/2) 
and therefore c R = c s = cos 6. With this and a coupling constant of 
k = mti/\ we obtain from equation (41) for the signal amplitude 

S = -j exp (-ad/cos fl-e^-sin (i> 2 + ft*/(l + $>')* 

v = irn^d/X cos 6 (46) 

£ = M/2 cos 6 = A0-/3dsin 0„ 

where v and £ have been redefined, and £ has been expressed in the needed 
form with the use of equations (14) and (18). 

Equation (46) has a form similar to that of equation (42) except 
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for an additional exponential term containing the absorption con- 
stant a. This term decreases the peak efficiency and it changes the 
angular sensitivity of the grating. But this change is very small, even 
for high loss values, as illustrated in Fig. 7. This figure compares the 
angular sensitivitities of a lossless grating (D = 0) with that of a 
grating of high loss (Z)„ = 2) for a parameter value of v = ir/2, a Bragg 
angle of O = 30°, and an optical grating thickness of (3d = 2wnd/X = 50. 
The loss parameter D was defined as 

D = ad/cos d (47) 

which is closely related to the conventional photographic density D 
(except that D is measured in the direction of the reference wave 
given by O ). A value of D = 2, which is the parameter used for the 
dashed curve, represents very high loss, with a decrease of the peak 
efficiency by a factor of about 50. Still, the differences of the two sen- 
sitivity curves are very small and consist mostly of an angular shift. 
The differences are even smaller for larger values of fid (we checked 
up to fid = 200), and, of course, for smaller values of D . The main 
conclusion is that the presence of loss has very little influence on the 
angular sensitivity of a dielectric transmission grating. This is probably 
because absorption influences the phase relations between the waves 
R and S very little. It agrees with observations by Belvaux. 31 

Next let us consider the influence of loss on the efficiency of a slanted 
dielectric grating. For simplicity we assume Bragg incidence, that is, 
# = 0. The obliquity factors are positive and given by c R — cos O 
and c s = —cos (0 O — 20). For this case we can write equation (41) 
for the signal amplitude S in the form 

S = -i^-exp [-*D (1 + c)] sin (u 2 - fl*/(l - €W 

V = TTtt^/XfcflCs) 1 (48) 

e = *D (i - c) 

where we have used the loss parameter D as above in equation (47), 
and the slant factor c 

D = ad/c n = ad/ cos 6 

c = Cr/c s = —cos o /cos (0 O — 2(f>). 

Figure 8 shows the diffraction efficiency of slanted grantings as 
calculated from equation (48). The efficiencies are plotted as a function 
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Fig. 8 — Transmission holograms — the efficiency of lossy dielectric gratings as a 
function of slant for v = 7r/2. c = c R /cs is the slant factor. 

of the slant factor c for various values of D , and for a value of v = ir/2 
which corresponds to the maximum attainable efficiencies. Similar 
curves for v = 7r/4 and v = St/4: and the same D values are almost 
identical to the curves of Fig. 8, except that the efficiency scale is 
reduced to a maximum efficiency of 0.5. This implies that for the 
range of chosen parameter values the exponential factor in equation 
(48) dominates in predicting the slant-dependence of the diffraction 
efficiency. 

The results show that, for higher efficiencies, the grating prefers 
small c-values, assuming constant 6 and D . This is a preference of 
small exist angles for S which means that we get the best efficiency 
if the signal wave leaves the grating on the shortest possible path 
after it has been generated. 

3.3 Unslanted Absorption Gratings 

When one records holograms in conventional photographic emul- 
sions one produces absorption gratings (bleaching can convert this 
into a dielectric grating). In an absorption grating there is no spatial 
modulation of the refractive index (n t = 0) and the coupling is provided 
by a modulation (ati) of the absorption constant. We have, then, an 
imaginary coupling constant k = —jai/2. In this section we study the 
efficiencies and the angular and wavelength sensitivities of unslanted 
absorption gratings where <j> = ir/2 and c R = c s = cos 0. From equa- 
tion (41) we obtain for the signal amplitude 



WAVES IN THICK HOLOGRAMS 2927 

8 = -cxp (-ad/c R ) -e- i( -sh (u 2 - f)»/d - WW 

v = ciid/2 cos 6 (49) 

£ = iM/2 cos « A0-/3d-sin O = -£(AX/X)2Cdtan 0„ 

where ^ and £ are real-valued, and equation (18) was used to express 
the parameter £ again in various forms, showing explicitly the angular 
deviations A0 and the wavelength deviation AX from the Bragg 
condition. 

For Bragg incidence we have £ = 0, and obtain from the above a 
formula for the diffraction efficiency 77 of absorption gratings 

77 = exp (-2ad/cos o )-sh 2 (M/2 cos 0„). (50) 

As we exclude the presence of negative absorption (gain) in the medium, 
there is an upper limit for the amplitude «i of the assumed sinusoidal 
modulation, which is a 1 ^ a. The highest diffraction efficiency possible 
for an absorption grating is reached in the limiting case a t = a for a 
value of ad/cos O = In 3. According to equation (50) this maximum 
efficiency has a value of i\ max = 1/27, or 3.7 percent. 

Figure 9 shows values for the diffracted amplitude S of absorption 
gratings as computed from equation (50) as a function of the modula- 
tion amplitude a t and for various values of the depth of modulation. 
For convenience we have again used loss parameters, which are D = 
ad/cos O and Dx = a,d/cos 6 . D x is a measure for the amplitude 
of the spatial modulation and D /Di = a/a x indicates the modulation 
depth. The dashed curves for constant D show the grating behavior 
for constant background absorption. We have plotted S on a linear 
scale in order to identify the regions of linear grating response. Note 
that a good linear response and relatively good efficiency is obtained 
if the absorption background is held constant to a value of about D = 1 • 

Equation (49) predicts also the angular sensitivity and the frequency 
sensitivity of absorption gratings. Such sensitivity curves are plotted 
in Fig. 10 for the special case of «i = a and values of v = Di/2 = 1 
(dashed) and v = \ In 3 = 0.55. For the latter parameter value the 
peak efficiency of 3.7 percent is reached, and for v = 1 we have a peak 
efficiency of 2.5 percent. In the figure the relative efficiencies are plotted 
as functions of the parameter £. We note that there is very little dif- 
ference between the sensitivity curves for the two v- values chosen. 
We have also computed the sensitivity for smaller values of v (0.2, 0.4), 
but the resulting curves differ so little from the ones shown that we 
have omitted them from the figure. The sensitivity curves are very 
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Fig. 9. — Transmission holograms — the diffracted amplitude of an absorption 
grating as a function of the modulation D t = aid/cos 6 = 2v for various modulation 
depths Di/Do (solid curve) and various bias levels Do = ad/cos (dashed curve). 



similar to those of the dielectric gratings with smaller ^-values which 
are shown in Fig. 6. Again, the half-power points are reached for about 
f = 1.5. But for absorption gratings there is no narrowing with in- 
creasing values of v, and the side lobe intensity remains low. 

3.4 Slanted Absorption Gratings 

In this section we consider the influence of slant on the efficiency 
of an absorption grating. For simplicity we assume Bragg incidence 
(t? = 0), and describe the slant by the obliquity factors c R = cos 9 
and c s = cos (0 O — 2(f)), as before. Using equation (41) we obtain, 
for this case, the following expression for the signal amplitude S 



S = - 



y'expf-i^i+yjsh^+rjva+fw 
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v = a,rf/2(c a Ca)* 



(51) 



where v and £ are redefined as real parameters. We have plotted the 
slant-dependence of absorption gratings in Fig. 11 for the special case 
of a, = a, that is, maximum depth of modulation. The diffraction 
efficiency -q is shown as a function of the slant factor c for various 
values of the loss parameter D . These quantities are defined, as 
before, by 



and 



D = ad/Cn = cad/ cos O 



c — c R /c s . 



(52) 



The efficiency is seen to reach its absolute maximum of 3.7 percent 
for the unslanted grating (c = 1) and for a loss parameter of D = In 3. 
For larger values of D the efficiencies reach relative maxima for exit 
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Fig. 10 — Transmission holograms — the angular and wavelength sensitivity of 
an absorption grating for a\ = a (Di = Do) and values of v = Di/2 = 0.55 (r) = 
0.037) and v = D,/2 = 1 (vo = 0.025). 
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Fig. 11- — Transmission holograms — the efficiency of an absorption grating as a 
function of slant for <*i = a (Di = Do), c = c R /c s is the slant factor. 



angles of the signal wave which are smaller than that of the reference 
wave (c < 1), while for smaller Z) -values the situation is reversed. 

3.5 Mixed Gratings 

Mixed gratings are those in which both the refractive index and 
the absorption constant are spatially modulated. This may occur in 
some recording materials (for example, as a result of incomplete bleach- 
ing, or in cases where strong absorption peaks are developed which 
cause refractive index changes according to the Kramers-Kronig rela- 
tions).* Mixed gratings are described by a complex coupling constant, 
which is given in equation (9). For the special case of unslanted gratings 
(<f> = 7r/2) and Bragg incidence (# = 0) equation (41) simplifies to 



S = —j exp (—ad/ cos O ) sin («d/cos O ) 



(53) 



where k is complex. From this we obtain, after some algebra, an expres- 
sion for the efficiency of mixed gratings 

r\ = SS* = [sin 2 (irriid/X cos O ) + sh 2 (at/2 cos O )] exp ( — 2ad/cos O ), 

(54) 

where n x and a x are the amplitudes of the modulation of the refractive 
index and the absorption constant, and a is the average absorption 



* Such effects have recently been observed by Nassenstein (see Ref. 32). 
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constant. We note that, at least for the special case considered here, 
there is a simple addition of the intensities diffracted by the dielectric 
grating and the absorption grating respectively [compare equations (46) 
and (50)!]. The exponential factor including a insures that the formula 
does not predict efficiencies larger than 1. 

IV. REFLECTION HOLOGRAMS 

In reflection holograms the recorded fringe-planes are of an orienta- 
tion which is more or less parallel to the surfaces of the recording 
medium, and the signal appears as a "reflection" of the reference wave. 
We have illustrated this situation in Fig. 4b. It is expressed in the 
coupled wave analysis by negative values of the obliquity factor 
Cs(cs < 0). In addition, the signal amplitude S of interest is obtained 
by evaluating the signal wave in the plane z = 0, which is also the 
entrance plane for the reference wave R. For reflection holograms a 
slant angle <f> = describes the case of unslanted gratings. Apart from 
these differences the following discussion of the detailed behavior of 
reflection holograms proceeds in a pattern similar to that of Section III, 
where we have discussed transmission holograms. 

From equations (30) and (39) we obtain a general formula for the 
signal amplitude of reflection holograms which can be written in the 
form 

8 = F*J -sh (y ch a)/ch (a + v ch a) 

v = JKd/\(c R c s ) h 

(55) 

, l (a a .&' 

\V R L S 0,5 

sh a = !/V 

where we have again defined (complex) parameters v, £ and a, which 
lump together the constants of the medium (n, a, n x ,a lf K) t the obliquity 
factors c R and c s , the wavelength, the grating thickness d and the 
dephasing measure &. 

4.1 Lossless Dielectric Gratings 

The lossless dielectric grating is characterized by a real-valued 
coupling constant k = irrii/X, and by zero absorption a = a, = 0. 
As in the transmission-hologram counterpart, it is easy to include the 
case of slant in the analysis. For the present case we can rewrite equa- 
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(56) 



tion (55) in the form 

8 = {$/{&> + (i - 6W-«* (/ - f) k \ 

v = JKnid/Xfpi&s)* 

£ - -M/2c a 

where v and £ have been redefined as real-valued parameters (c s is 
negative!). 

The associated formula for the diffraction efficiency of lossless di- 
electric gratings is 

v = i/{i + (l - *W>M" O' - tV). (57) 

which also provides a description of the angular and wavelength sen- 
sensitivities of the grating. For unslanted acoustic gratings this formula 
has been previously given by Quate and his associates. 22 Sensitivity 
curves calculated from equation (57) are shown in Fig. 12, where the 
normalized efficiencies are plotted as a function of £ for various values 
of v = const. The figure shows the sensitivity of a grating with v = tt/4 
and a peak efficiency of 43 percent, a grating with v = ir/2 and i) = 0.84, 




Fig. 12 — Reflection holograms — the angular and wavelength sensitivity of a 
lossless dielectric grating with the normalized efficiency 77/770 as a function of £. 
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and the corresponding values for v = 3ir/4 and t} = 0.96. For v = 7r/4 
the half-power points of the grating response are reached for values of 
approximately £ = 1.7. But there is considerable broadening of the 
sensitivity curves for increasing values of v, and an increase in the 
side-lobe level. 

As in equation (44) for transmission holograms, we can express the 
parameter £ directly in the angular deviation A0 or the wavelength 
deviation AX by using equation (18) to obtain 

£ = Ad-Kd-s'm (0 O - 4>)/2cs 

= A\-K 2 d/8wnc s . (58) 

These expressions can again be used to formulate rules for the angular 
bandwidth and the spectral bandwidth of the grating. 

For an unslanted grating (<£ = 0) and Bragg incidence we have 
Cr = — Cs = cos d , and equation (57) simplifies to 

7] = th 2 (TTWxd/X cos do). (59) 

This is a formula which has been obtained previously for light diffraction 
by acoustic waves. 33,22 

4.2 Lossy Dielectric Gratings 

Let us first discuss the influence of loss on the angular and wave- 
length sensitivity of unslanted dielectric gratings. Here we have <t> = 
and, to a good approximation 

c R = cos O (1 — Ad tan O ) = cos 6 

c s = -cos 0„(1 + A0tan O ), (60) 

= - cos 0(1 + 2AX/X) 

at least as long tan 0„ ^ 1. One can show that the formula for the signal 
amplitude S, which we have given in equation (56), is still applicable 
for the present case of an unslanted lossy grating if we modify the 
parameters v and £ to 

v = irriid/X cos O 

£ = £ - jD , (61) 

So = -A0-/3dsin O 

D = ad /cos O 

where £ is now a complex parameter with £o describing the angular 
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deviations and D representing the loss. An evaluation of this formula 
is shown in Fig. 13, which shows the angular sensitivity of dielectric 
gratings for various values of the loss parameter D and a grating 
parameter of v — t/2. In constrast to what we have observed in the 
case of dielectric transmission holograms (Fig. 7), we see here a quite 
noticeable effect of the grating loss on the sensitivity curves. With 
increasing loss values the curves broaden in the wings, sharpen some- 
what in the center and the side-lobe level decreases. 

To study the influence of loss on the diffraction efficiencies of dielec- 
tric gratings we rewrite equation (55) in the form 



S = i^j/Wv + (l + OV-ooth (p 2 + ft*] 



v = JTrnid/XfpuPs)* 
£ = i£> (l - c) 



(62) 



where we have written v and £ as real-valued parameters in a form 
which is valid for Bragg incidence and for slanted or unslanted gratings. 
Just as in the case of transmission holograms we have used the loss 




Fig. 13 — Reflection holograms — the influence of loss on the angular and wave- 
length sensitivity of a dielectric grating for v = ir/2. The normalized efficiencies 
v/vo are shown. The peak efficiencies are jjo = 0.84 for Do = 0, vo = 0.64 for 
Do = 0.5, jfo = 0,28 for D = 1, and 70 = 0,12 for D = 2, 
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parameter D„ and the slant factor c (which is now negative) 

D = ad I COS O (nn\ 

c = c R /c s . 
In the case of unslanted gratings the parameters v and £ simplify to 



v = 7rn,rf/X cos 6 

£ = D = ad/cos O • 



(G4) 



The results of a numerical evaluation for unslated gratings are shown 
in Fig. 14, where the signal amplitude is plotted as a function of v for 
various values of the loss parameter D . The curve D = gives the 
values for lossless gratings, while the others indicate the influence 

of loss. 

The behavior of slanted dielectric gratings in the presence of loss 
is shown in Fig. 15. The curves of this figure are also computed from 
equation (62) and show the diffraction efficiency as a function of the 
slant factor for v = tt/2 and various values of the loss parameter D . 
For constant D we notice an increase of the efficiency for decreasing 
values of the slant factor, as in the case of transmission holograms. 




0.4 0.6 0.6 

v/ir - n,d/\cos0 o 



Fig. 14 — Reflection holograms— the influence of loss on the diffracted amplitude 
S of an unslanted dielectric grating. | S | is shown as a function of v/w = nid/\ cos 6 
for various loss parameters Do. 
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Fig. 15 — Reflection holograms — the efficiency of a lossy dielectric grating as a 
function of slant for v = v/2. c = c R /c s is the slant factor. 



Again, for given loss and a given angle of incidence short signal paths 
through the grating (that is, small exit angles) are preferred for higher 
efficiencies. 

4.3 Unslanted Absorption Gratings 

Following the pattern set in the discussion of transmission holograms 
(Section III), we again describe an absorption grating by an imaginary 
coupling constant k = — ja : /2, and proceed to study the diffraction 
efficiencies and the angular and wavelength sensitivities of unslanted 
(<t> = 0) gratings. In this case equation (55) simplifies to 

s = -;(jr)V{*A + W»? - i]*coth & - o»| 

v = ja^McnCs? (65) 

£ = -Do - j£o 

where the real-valued parameters D and £ can be expressed to first 
order in the angular deviations A 6 and the wave-length deviations 
AX by 

Do = ad/ cos d . 

So = Ad-0ds'm d Q = i(A\/\)Kd. 
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D is a loss parameter as before, and £ is a normalized measure for 
the angular or the wavelength deviations from the Bragg condition. 

If the Bragg condition is obeyed equation (65) can be written in 
the form 

S = -D,/2[D + (D n 2 - D?/4)*-coth (Dl - Dj/4) 1 ] (67) 

where 

Di — 2v x = aid/cos O 

measures the spatial modulation of the absorption constant (c*i). 

For the deepest allowable modulation where we have D y = D (ai = a ), 
this equation predicts the maximum diffraction efficiency 77 mnx which 
is possible for reflection holograms with a (sinusoidal) absorption 
modulation. We obtain 77 mBX = 1/(2 -f- V3) 2 , or a maximum efficiency 
of 7.2 percent for D = Di—> °° . The formula reflects the experimental 
fact that, for reflection holograms of the absorptive kind, one obtains 
the largest efficiencies for high photographic densities. Figure 16 shows 
a numerical evaluation of the above formula. Here the signal amplitude 
S is plotted as a function of the modulation amplitude D t for various 
levels of loss "bias" D (dashed curves) and for various modulation 
depths Do/Di . 

An evaluation of the grating sensitivity as predicted by equation (65) 
is shown in Fig. 17 for the special case of a maximum depth of modula- 
tion where Z), = D . In this figure the (normalized) efficiency is plotted 
as a function of the parameter £ for various values of Z) x = D . As in 
the corresponding grating for the case of transmission holograms (Fig. 
10) the sensitivity curves are seen to reach their half -power points for 
values of about £ = 1-5- But in the present case of reflection holograms 
there is a noticeable broadening of the curves with increasing loss 
values D, = D . 

4.4 Slanted Absorption Gratings 

In this section we consider the influence of slant on the diffraction 
efficiency of an absorption grating for reflection holograms. We assume 
Bragg incidence (t? = 0) and again use the obliquity factors c R = cos 0„ 
and c s = — cos (0„ — 2<f>) to describe the slant (for reflection holograms 
we have c. s < !) . We find that equation (65) can be used as a formula 
for the signal amplitude for the present case if we modify the parameters 
to 

v = ;« 1 rf/2(c R c 5 ) i = | /)»(<:)* 
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Fig. 16 — Reflection holograms — the diffracted amplitude of an absorption 
ating as a function of the modulation A = aid/cos = 2v for modulation depths 
/Do (solid curve) and bias levels Do = ad/cos 6 (dashed curve). 



f = *D (1 - c) 

D = ad/ cos O , A = «id/cos O 

C = C ff /Cs 



(68) 



where the slant factor c is negative. All these parameters are real-valued 
in the present case. For a maximum depth of modulation, that is, a x = a, 
there are further simplifications, and we obtain a simple expression for 
the slant-dependence of the diffraction efficiency 

,, = - C /{1 - c + (1 - c + c 2 ) J -coth i/> (l - c + c 2 ) 1 } 2 . (69) 

Figure 18 shows a numerical evaluation of this formula for various 
values of D Q = D Y . The slant factor value of | c | = 1 refers to unslanted 
gratings. In this case the maximum efficiency value rj max = 0.072 is 
approached for large D . We note that for values of D„ below unity 
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the efficiencies increase for | c | -values larger than 1 and up to about 3, 
that is, for relatively large exit angles of the signal wave. 

4.5 Mixed Gratings 

Mixed gratings are described by a complex coupling constant 
K = 7rn,/X — ;ai/2 (see Section 3.5). For Bragg incidence (# = 0) 
and unslanted fringe-planes (<£ = 0) we can obtain from equation (55) 
a formula for the signal amplitude of mixed gratings, which is 



8 --*/{« + ^ + ^- coth c^^ + « 2)i } 



(70) 



where k is of complex value, a is the average absorption constant, d the 
grating thickness and O the angle of incidence. 

V. AMPLITUDES OF THE DIRECT WAVES 

For diagnostic purposes it is often of interest to monitor the change 
in amplitude of the direct reference wave R, which is depleted because of 
diffraction into S and absorption. The quantities of interest are the 
amplitudes R(d) which can be obtained from the analysis of Section 2.2. 




Fig. 17 — Reflection holograms— the angular and wavelength sensitivity of an 
absorption grating for «i = a (D, = D ) and values of Di = 2v = Do = 0.2 
(1,0 = 0.007), D, = 1 (no = 0.05), and D, = 2 (^ = 0.068). 
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Fig. 18 — Reflection holograms — the efficiency of an absorption grating as a 
function of slant for a t = a (Di = Do), c = c R /c s is the slant factor. 

We will give here the general results for transmission and reflection 
holograms. The notation is that of Section 2. 

5.1 Transmission Holograms 

From equations (27) and (33) we get for the constants r, of equation 
(25) the expressions 

Ti = -k 2 /Cs(7i — T2)(Cft7i + «) / 71 n 

r 2 = k/c s (Ti - 72)(c fi 72 + «)• 

Using this we can write the output amplitude R(d) of the reference 
wave in the form 

m = -r^—, i' s3Lh r 1 - e - m ¥ 1 ) ■ <") 

Cs (Ti — 72) \Cr72 + a C«7, + a/ 

5.2 Reflection Holograms 

For reflection holograms we use equations (27), (37), and (39) to 
express the constants r, in the form 

r, = (csTi + « + fi) exp (y 2 d)/{exp (y a d)(a + 38 + c s7l ) 

- exp (7id)(o + j& + c S 7 2 )} « 3 * 
r a = - (c 5 7 2 + a + #) exp (7,d)/{exp (73d) (a + # + c 5 7i) 

- exp (7,d)(a + it? + C S 7 2 )}. 
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The output amplitude R(d) of the reference wave becomes 
R(d) = c s (t, - 72)/ {(« + ft + c S 7,) exp (-7^) 

- (a + ^ + C S 7 2 ) exp (-Tad)} • (74) 

More detailed evaluations of the above formulas should follow the 
pattern prescribed in Sections III and IV. They can be undertaken 
for the specific case when the need arises. 

VI. VALIDITY OF THE THEORY 

We have tried to make our results as generally applicable as possible. 
We have allowed for the presence of absorption in the various hologram 
gratings and for a slant of the fringe planes. But a whole range of assump- 
tions had to be made to make the simple coupled wave analysis possible. 
It seems appropriate to recount these assumptions to make clear the 
region of validity of the coupled wave theory. We have assumed that: 

(i) The electric field of the light is polarized perpendicular to the 
plane of incidence. However, the appendix gives an extension of the 
theory to allow also for light of parallel polarization. 

(ii) A slant of the fringe planes with respect to the z-axis is allowed, 
except that these planes are perpendicular to the plane of incidence. 
(This is reflected in the assumption e(x, z), <t(x, z).) But this assumption 
is not made in the generalization which we have given in the appendix. 

(Hi) The spatial modulation of the refractive index and the absorption 
constant is sinusoidal. 

(iv) There is a small absorption loss per wavelength and a slow 
energy interchange (per wavelength) between the two coupled waves. 
This condition is stated mathematically in equation (7) and justifies 
neglecting the second derivatives R" and S" in the analysis. 

(v) There is the same average refractive index n for the regions 
inside and outside the grating boundaries. If the grating has interfaces 
with air, then Snell's law has to be used to correct for the angular 
changes resulting from refraction. 

(pi) Light incidence is at or near the Bragg angle and only the diffrac- 
tion orders which obey the Bragg condition at least approximately are 
retained in the analysis. The other diffraction orders are neglected. 

A detailed mathematical justification of assumption vi is outside the 
scope of our simple analysis. One can advance physical arguments to 
show that this step limits the validity of the theory to "thick" gratings, 
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where the phase synchronism between the two coupled waves has 
enough time to develop a strong and dominating effect. Better definitions 
of a "thick grating" must come from more accurate theories which are 
available for special cases. A large amount of work has been done on 
acoustic diffraction gratings which correspond to the case of our 
unslanted, lossless, dielectric transmission-hologram gratings. 30 In 
acoustic diffraction one defines the parameter 

Q = 2ir\d/nA 2 (75) 

as an appropriate measure of grating thickness. We can regard a grating 
as thick when the condition Q ^> 1 holds. 21 ' 10 It appears that the coupled 
wave theory begins to give good results for values of Q = 10. This is 
particularly well demonstrated by Klein and his associates in theoretical 
and experimental work on acoustic gratings for the predictions of both 
the peak efficiencies and the angular sensitivities. 16,17,34 We hasten to 
add that for the majority of practical holograms the parameter Q is 
larger, and sometimes much larger, than 10. 

Further checks of the validity of the coupled wave theory are provided 
by comparisons with accurate computer calculations and with experi- 
ments on special examples of gratings. Burckhardt has made computer 
calculations on unslanted, lossless, dielectric transmission holograms 
for selected values of grating parameters which are commonly encoun- 
tered in holography. 18,19 Comparison with the results of the coupled 
wave theory shows very satisfactory agreement. 35 Measurements by 
Shankoff and Lin on dielectric transmission holograms prepared with 
dichromated gelatin yielded diffraction efficiencies approaching 100 
percent, which agrees with the theory (even though there may be some 
uncertainty as to the exact nature of the refractive index variations). 1 ' 

Efficiency measurements on thick absorption gratings for the case of 
transmission holograms were made by George, Mathews, and Latta. 36,37 
Efficiencies approaching our predicted maximum value of 3.7 percent 
were observed. 

Kiemle has studied unslanted (# = 0) reflection holograms for the 
special case of normal incidence (0 O = 0) by analyzing equivalent 
four- terminal networks. 38 His treatment of absorption gratings cor- 
responds to the material we discussed in Section 4.3 specialized to the 
case of O = 0. But Kiemle's value of 2.8 percent for the maximum 
diffraction efficiency of absorptive reflection holograms does not agree 
with our prediction of 7.2 percent. This disagreement appears to derive 
from a set of restrictive assumptions made in Kiemle's work. Experi- 
mental observations on absorptive reflection holograms were made by 
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Lin and Lo Bianco. 9 Efficiency values as high as 3.8 percent were 
measured, which seems to support the predictions of the coupled wave 
theory. But further experiments are needed for a good confirmation. 

VII. CONCLUSIONS 

We have discussed a coupled-wave analysis of the Bragg diffraction 
of light by thick hologram gratings. This approach made it possible to 
derive simple algebraic formulas for the behavior of various types of 
holograms, even for the case of high diffraction efficiencies where the 
incident wave is strongly depleted. The treatment covers transmission 
holograms and reflection holograms, and it includes the spatial modula- 
tions of both the refractive index and the absorption constant. The 
influence of loss in the grating and of slanted fringes is also discussed. 
Formulas and their numerical evaluations are given for the diffraction 
efficiencies and the angular and wavelength sensitivities of various 
grating types. 

For special cases we can compare the results of this theory with more 
accurate computations and with experimental observations. These 
comparisons give us the confidence to assume that the coupled wave 
predictions are good for a broad range of practical hologram types. 
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APPENDIX 

Reduced Coupling for Light Polarized in the Plane of Incidence 

In the body of this paper it was assumed that the incident light is 
polarized perpendicular to the plane of incidence. The purpose of this 
appendix is to show that we can use the results of the main paper also 
when the light is polarized in the plane of incidence, provided that we 
modify the coupling constant k. Such a modification is suggested already 
by the dynamical theory of X-ray diffraction. 

As in Section II we start with the wave equation 

V'E - V(V-E) + fe 2 E = (76) 

for the electric field in the grating. Here, in contrast to equation (1), 
we have described the field by the vector quantity E and have included 
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the term V(V-E), which is not necessarily zero. The constant k 2 is 
denned in equation (4). As in the main paper, we assume that only two 
waves are present in the grating, and put 

E = R(z)<f "* + S(e)e-'"* (77) 

using the vectors R and S to describe the amplitudes of the reference 
and signal waves, p and d are the propagation vectors (as in Section II) 
which point in the direction of the wavenormals. They are related by 
equation (11). In addition we assume that, both, R and S are transverse 
waves, that is, that the following conditions hold 

(e-R) - o, (78) 

(d-S) m 0. 

Combining equations (76), (77), and (78) we get, after separating 
terms with equal exponentials and neglecting second derivatives d 2 /dz 2 

-2j Pz R' + jgR: - 2j«/3R + 2k/3S = (79) 

-2j<r,S' + jdS', + (/3 2 - <r 2 - 2ja0)S + 2k/3R = (80) 

where R. and S M are the 2-components of R and S, and the notation of 
Section II is used. 

We now make the additional assumption that the polarizations of 
R and S do not change in the grating and write 

R(*) = R(z)i, (81) 

S(z) = S(z)s, 

where R(z) and S(z) are the scalar amplitudes of the two waves, and 
r and s are polarization vectors independent of z. These vectors are 
normalized so that 

(r-r) - 1, (s-s) = 1. (82) 

Because of (78) we have 

(r- e ) = 0, (sd) = 0. (83) 

After forming the dot products of r with eq. (79) and of s with (80) 
we use equations (81), (82), and (83) to arrive at 

-2j Pt R' - 2ja(3R + 2k0S(t-s) = (84) 

-2j<j,S' + (P 2 - a 2 - 2jafi)8 + 2k0R(t-s) = 0. (85) 
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As in Section II, we introduce the abbreviations 

c„ = p./0, c, = a./fi, (86) 

and 

d = (f - a 2 )/20, (87) 

which allow us to write the above equations in the form 

c R R' +aR = -jk(t-s)S (88) 

c s S' + (a + j&)S = -jk(t-s)R. (89) 

These are coupled wave equations which govern the Bragg diffraction 
of light polarized parallel to the plane of incidence, and indeed, of light 
of arbitrary polarization. They are similar in form to the coupled wave 
equations (21) and (22) which were derived for perpendicular polariza- 
tion. The only difference is a reduction of the effective coupling constant 
by the dot product (r • s) of the two polarization vectors. 

Referring to the grating geometry of Fig. 1 we have (r-s) = 1 for 
light polarized perpendicular to the plane of incidence. For parallel 
polarization the value of this dot-product depends on the inclination 
angles, and we have a reduced effective coupling constant k b given by 

K , = »c(r-s) = -K COS 2(«o - 4)- ( 9 °) 

We can apply the results of the main paper for parallel polarization if 
we replace k by k, . For this polarization there is the trivial case of a 
Bragg angle of 45° (that is, diffraction angles of 90°) where (r-s) = 
and the intensity of the diffracted light goes to zero. 
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